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Abstract 

in 
o 

■ We derive the zero mode solutions for a Majorana fermion in the background of 

a cosmic string and contrast it with the zero mode solution for a (neutral) Dirac 

Ch ! 

fermion. A Majorana zero mode carries no vector or axial charge, and it cannot 

^_ 

be bosonised. We study the implications for vorton formation and stability. In 
the massless limit stability of the zero mode is guaranteed by energy-momentum 
^NJ ■ conservation. However, zero modes obtain an effective mass on string loops. It is 

| found that the conditions under which current formation can be effective are exactly 

in ' 



those for which zero mode decay is most likely to occur. 



1 Introduction 



The existence of Majorana mass terms is particularly important for particle-physics phe- 



nomenology in the context of neutrino physics [I]. Indeed, the see-saw mechanism j2], 
which can naturally explain the smallness of the observed left-handed neutrino masses [3. , 
relies on the existence of a Majorana mass term for a standard model singlet, the right- 
handed neutrino. Right-handed neutrinos and the see-saw mechanism are a prediction 
of unified theories which are left-right symmetric, such as Patti-Salam and SO (10) grand 
unified theories (GUTs). After breaking of the gauged U(1)b-l symmetry a Majorana 
mass for the the right-handed neutrino is generated through the Higgs mechanism. 

In the context of cosmology, GUTs predict the existence of cosmic strings jUIHj- If the 
Higgs field which form the string breaks U(1)b-l, B — L cosmic strings are produced [HllZj- 
They can be important for baryogenesis |H]. Cosmic strings can also contribute to 
the cosmic microwave background anisotropies and they produce gravitational waves 0. 
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Further, fermions coupling to the string forming Higgs field can be trapped on the string, 
travelling along the string at the speed of light ^0] • Such fermionic zero modes can give 
rise to superconducting strings jllj . Conserved fermionic currents can stabilise cosmic 
string loops against gravitational collapse, resulting in the formation of vortons which 
are potentially cosmologically catastrophic It was shown, however, that if there are 
fermionic zero modes which travel in opposite directions on the string they could scatter 
with each other, thereby strongly reducing the current |T?H I14j. Hence, string loops with 
chiral currents are the best vorton candidates. In an anomaly free theory, chiral zero 
modes have to be neutral; they can either be Dirac or Majorana. 

In this paper we will be interested in chiral zero modes and in particular in zero 
modes for Majorana fermions. Such zero modes arise naturally in B — L cosmic strings, 
as the U(1)b-l breaking Higgs field gives a Majorana mass to the right-handed neutrino 
They can produce the observed lepton asymmetry of the universe, independently of 
the reheating temperature after inflation [HI U\ ■ We will discuss the difference between 
Majorana and Dirac zero modes. Neutral Dirac currents are persistent by virtue of a 
conserved quantum number, whereas Majorana currents are persistent for kinematical 
reasons. The difference between the two currents shows up when the zero-mode obtains an 
effective mass. Decay and annihilation of Dirac zero modes can be suppressed or forbidden 
by virtue of lepton number conservation; no such suppression exists for Majorana fermions. 
We apply our results to vortons. Neutrinos zero modes were discussed before in a slightly 
different context |15j . 

This paper is organised as follows. In the next section we will give a brief introduction 
to Majorana and Dirac spinors. In section El we present the zero mode solutions for both 
the Dirac and Majorana spinors in the presence of a Nielsen-Olesen string. We quantise 
the string in section and compute the conserved currents. Moreover, a description of 
the effective theory in 2D is given. We apply our results to vortons in section |3J Both 
current condensation and stability is discussed, highlighting the difference between Dirac 
and Majorana currents. We end with some concluding remarks. 

2 Dirac and Majorana spinors 

In the early Universe when particles are still massless, there is no distinction between 
Dirac and Majorana fermions. All fermions and anti-fermions can be described by two- 

2 



component Weyl spinors. In grand unified theories fermions and anti-fermions are in the 
same representation and therefore must be described by spinors with the same chirality. 
By convention, we use left-handed Weyl spinors. 1 Weyl spinors in 4-dimensions have two 
on-shell degrees of freedom. 

Fermions get masses via the Higgs mechanism. The mass term involves a coupling 
between a fermion bilinear and a scalar field <fi which acquires a non-vanishing vacuum 
expectation value (VEV). The mass term must be invariant under all symmetries of 
nature, and in particular, it must be invariant under Lorentz transformations. 

Given two left-handed Weyl spinors ipL and xl, the Lorentz invariant bilinear which 
can be constructed is of the form 

xIc^l^ylcTl ■ (1) 

If the Weyl spinors transform under some symmetry, the above term is in addition gauge 
invariant if ipL and xl are charge conjugate of each other, i.e., if Xl — "0£- The bilinear 
term given in Eq. then becomes: 

m T C*pL + h.c. = (CTrYC^l + h.c. = XbP T C X l + h.c. = Tr^l + Mr- (2) 

This leads to the standard Dirac mass term 

£d = -if)Rm D ip L - ipLin^ipR (3) 

T 

with rrif) = A0 and A a Yukawa coupling. The two Weyl spinors ipi and i/jr = Cipl 
combine to make a Dirac spinor 

^Pd = ^l + *Pr, (4) 

and if ttid is real, the mass term can be written as Lb = —m pip Dip £>. A Dirac spinor 
has on shell four real degrees of freedom corresponding to particle/antiparticle and pos- 
itive/negative helicities. The Dirac mass term is invariant under a global U(l) transfor- 
mation; this implies the existence of a quantum number (such as lepton number) which 
enables to distinguish particles from antiparticles. A Dirac mass term can be written for 
both neutral and charged particles. 

1 Our notation and conventions are detailed in the appendix. 
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Now starting with a single neutral Weyl spinor, i.e., a particle which does not carry 
any quantum number under any unbroken symmetry, we can also construct, by setting 
Xl = ?Pl m Eq. (Q), the following Lorentz invariant: 

^IC^L + TlCTl = Wk^L + Titfh. (5) 
We can then write down what is called a Majorana mass term 

£m = -^ c R m M ^L ~ ^Lm M ^ c R , (6) 

where tum — A0 and the factor half is conventional. Analogously to the Dirac spinor, a 
4-component Majorana spinor can then be defined 

IpM = IpL + Clp c L T = l/j L +ljj R = l/j M ,L + 1pM,R- (7) 

If itim is real, the mass term simplifies to Cm — ~\ m M'4 ) M'4 , M- A Majorana mass term is 
not invariant under a global U(l) transformation, and there is no quantum number which 
allows to distinguish particles from antiparticles (for instance, lepton number cannot be 
defined for a Majorana neutrino). A Majorana spinor is made of a single Weyl spinor 
and it has on shell two real degrees of freedom corresponding to the two helicity states. 
It satisfies the Majorana condition 

1>M = ^M, (8) 

that is, a Majorana particle is its own antiparticle. 

Majorana masses are particularly important in understanding the smallness of the 
neutrino masses via the see-saw mechanism Indeed, the latter relies on the existence 
of a superheavy Majorana mass term for the right-handed neutrino on top of a standard 
Dirac mass term which arises from the coupling of both left and right-handed neutrinos 
to the standard model Higgs. In principle there can also be a small Majorana mass 
for the left-handed neutrino. There are then two distinct mass eigenstates which are 
both Majorana: a state with a tiny mass which is a superposition of mostly left-handed 
neutrino with a small admixture of the right-handed neutrino, and a superheavy state 
which is mostly composed of the right-handed neutrino. In the context of topological 
defects, we will be interested in a Majorana mass for the right-handed neutrino which 
arises as a consequence of gauge symmetry breaking. 
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Table 1: Charges of the Higgs field and of the fermion fields ipL and xl, which acquire 
respectively Dirac and Majorana mass by coupling with 0, under the broken U(l)x and 
any unbroken U(1)q. Here r, n ^ 0, but q can either be zero or non-zero. 



In summary, massive charged fermions are always Dirac fermions whereas neutral 
fermions can either be Dirac or Majorana. Which mass terms can be constructed depends 
on the quantum numbers of the fermions and the scalar fields of the model. To be 
explicit, consider a U(l)x x U(1)q gauge symmetry, possibly embedded in a larger gauge 
group. The U(l)x symmetry is broken by the non-vanishing vacuum expectation value 
of a complex Higgs field 0, while U(1)q remains unbroken. Given two left-handed Weyl 
fermions ipi and xl, Dirac and Majorana mass terms involving the bilinear term given 
in Eq. (0) multiplied by the Higgs field can be constructed if the charges of the fields 
are as given in Table 1. Let's call X and Q the charges of the fields under U(l)x and 
U(1)q respectively. breaks U{l)x and keeps U{\)q unbroken, hence X(0) = r ^ and 
Q(4>) = 0. Let X(iPl) = n and X(xl) = P- Then the Yukawa term is gauge invariant 
only ifn + p + r = and Q(?Pl) = —Q(xl)- We distinguish three cases: 

1. If Q(?Pl) = —Q(xl) 7^ 0, we have a Dirac mass term. The spinors ipL and xl 
combine to form a Dirac spinor which has four on shell degrees of freedom. It describes 
a charged fermion. 

2. If Q(ipi) = and n = —r/2, we have a Majorana mass. A Majorana spinor is 
formed out of x/jl, it has two on-shell degrees of freedom. It describes a Majorana fermion, 
a particle which is its own antiparticle. 

3. Now if Q(ipi) = and n ^ r/2, we have a Dirac mass term for a neutral particle. 
In this case, ipL and xl combine to form a Dirac spinor. Note that in this case always a 
global charge can be defined which distinguishes particles from antiparticles. If there is 
a second Higgs field with charge X = —2X(ip), respectively X = — 2X(x), a Majorana 
mass term can be constructed for ip (respectively x)- There are then two distinct mass 
eigenstates which are Majorana. 
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3 Chiral currents 



3.1 Zero modes solution 

In this section, we write down the fermionic equations of motion in the background of a 
cosmic string. We give the zero mode solutions for both Dirac and Majorana fermions, 
and compare the results. We only discuss Abelian cosmic strings; our results can be 
generalised straightforwardly to non- Abelian strings. 

Consider a U(l)x gauge symmetry which is broken down to the identity by a Higgs 
field 0, when the latter acquires a non-vanishing VEV. Since TTx(U(l)x) 7^ 0, cosmic 
strings form at the U(l)x breaking scale according to the Kibble mechanism [3]. For a 
straight infinite cosmic string lying along the z-axis, the Higgs field and gauge field A 
have the form 

= M\r)e in \ (9) 

A 6 = -n^-, A z = A r = 0, (10) 

er 

with 0o the Higgs VEV in the vacuum, and n G N the string's winding number. In 
the following we will concentrate on the minimum energy configuration corresponding 
to \n\ = 1. The profile functions f(r) and g(r) must satisfy the following boundary 
conditions 

/(0) = <7(0)=0, (11) 
lim f(r) = lim g(r) = 1. (12) 

r— >oo r -^oo 

The exact form of the functions f(r) and g(r) depends on the Higgs potential 0, which 
we do not specify. 

We first discuss the well known zero mode solution ^U] for a Dirac fermion getting its 
mass by coupling to 0* [TT] . The Lagrangian is 

C D = Jli^D^L + Jd^D^R - A0*^Vl - A0V^0*, (13) 

with = dfj, + iqAfj,, X the Yukawa coupling, and q the charges of ipL,R under U(l)x- 
Note that q<f> — qi + cLr = 0; in general qi and — q^ can be different. The equations of 
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motions are: 



(14) 



From Jackiw and Rossi ^U], and an index theorem we know that there is one nor- 
malisable zero mode solution ip° (a static solution) which is z and t independent. The 
zero mode solution are eigenvectors of the projection operator ^T] 



(15) 



with n = 1 for a vortex and n = — 1 for an anti- vortex. With our choice of gamma 
matrices (see the appendix) the zero mode solutions are 





a+(r,i 




\ 



/ 



(<x-(r,6)\ 



V/3-(r,C)y 



(16) 



for a vortex and an anti- vortex respectively. The wave functions a±(r, 0) and /3±(r, 0) are 
of the form f(r)e tW with 1 an integer and f(r) a real function localised on the string. We 
define 



0± 



rdrd# a±(r, 9), 
rdrd6f3±(r,6). 



(17) 



The ratio a±/ j3± is determined by the equations of motion, more specifically, by the 
charges qi and q^. Further the sum {a± + 0±) is set by the spinor normalisation. Hence 
there is no freedom in choosing the (relative) magnitudes of a± and (3±. We normalise 



The 4-dimensional solutions of the Dirac equation are of the form 2 



i) D = ip n e 



±iE(t-nz) 



(18) 



(19) 

the zero mode 



2 For a Dirac particle, let's call it rj, getting its mass by coupling to <j> instead of ■ 
solutions are f?„ = j = i>n=—i and rj^ l= _ 1 — ipn=i- It follows that for a vortex a fermion zero mode due to 
Yukawa coupling to <p moves in opposite direction as a fermion zero mode due to Yukawa coupling to (jf . 
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where we have set E = k = k z , with k > 0, for a string aligned along the z-axis. The 
zero-modes are massless excitations which travel at the speed of light along the string. 
Fermions trapped on a vortex (anti-vortex) travel in the +z (—z) direction. Note that 
both particles and anti-particles move in the same direction along the string, the direction 
being determined by whether the Yukawa coupling in Eq. (j!3|) is to a vortex or anti- vortex. 
The projection operator in Eq. (|15jl singles out one helicity state (determined in terms 
of a + and f3 + or a~ and /?_). The zero mode solution has thus only two real degrees of 
freedom which correspond to particle and anti-particle. It is effectively a 2-dimensional 
Weyl fermion. 

Now consider a Majorana fermion getting its mass by coupling to 0*. The Lagrangian 
reads 

C M = Tl^D^l - ~A<^|V>£ - ^WTl^r (20) 
The equations of motions are 

i^D^r = A0*Vl- (21) 

Normalising the U{\)x charge of (f> to unity, the charges of ipi and ip c R are qi = — = 
— 1/2. The two equations of motion are not independent, rather the second equation in 
Eq. (|21[) can be obtained from the first one by charge conjugation. The solution to the 
equation of motion is 

^™ = -=(V£(r, e)e- lE{t ~ nz) + tfj° n *{r, e) e +iE ^- nz )) (22) 
V 2 

with tp® as given in Eq. (fT^j) with the additional constraint that 

a± (r,9) = T P±(r,9). (23) 

The Majorana solution is just a linear combination of the Dirac solutions given in Eq. ()19j) 
which satisfies the reality constraint of Eq. (jHJ) . There is only one solution corresponding 
to one real degree of freedom, which is half the number of degrees of freedom of the 
Dirac solution. The zero mode is effectively a 2 dimensional Majorana- Weyl fermion A 
consequence of Eq. (123)1 is that both the vector and axial currents are zero. Nevertheless 
there is energy-momentum flowing along the string. 

8 



3.2 Quantisation and conserved currents 

In this section we quantise the Dirac and Majorana zero mode fields, and evaluate the 
conserved charges and energy-momentum. 

As before, we start with the Dirac spinor. We use the classical wave solutions to the 
Dirac equation, Eqs. (|TH| H9 j) . to write the quantised field 

{a(nk)r n e- ikit - nz) + b\nk)^J k ^) , (24) 

with as before k = k z . Note that k > always. Thus a(k < 0) are annihilation operators 
for states moving in the —z direction, which can only be exited for an anti- vortex (n = — 1). 
Likewise a(k > 0) are annihilation operators for states moving in the +z direction, which 
can only be exited for a vortex (n = 1). The wave function ip° is fc-independent by virtue 
of separation of variables in the equation of motion; the normalisation is also momentum 
independent. 

Imposing the equal time anti-commutation relations 

{ip(x, y, z),ift(x, y, z')} = 5{z - z^ipo^t (25) 

it follows that the operators a, b satisfy the anti-commutation relations {a(k), a^(k')} = 
{b(k), tf(k')} = 27i8(k — k'). One can construct a Fock space in the usual way [T7j . 

The vector current j M = ipj^ip and axial current j^ 5 = ip'j^'j 5 '^ are identically zero for 
H — 1, 2. The non-zero components for fi = 0, 3 are related, j 3 = nj° and j 35 = raj 05 , by 
virtue of Eq. (fTK|) . The conserved charges corresponding to the vector and axial currents, 
Q = f d 3 x( ) and Q 5 = f d 3 x( : j 05 : ) respectively, are 

Q = / 

Jo 
= N- 

Q 5 = / 

Jo 

= (a 2 

where, to avoid notational cluttering, we have dropped the subscript n to indicate whether 
it concerns a vortex or anti- vortex. Here : : denotes normal ordering, i.e., putting all 
annihilation operators to the left; this procedure automatically subtracts the zero point 
contributions. 
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^(aUk)a(k)-bUk)b(k)) 
- N, (26) 

° c\k 

— ({-a 2 + (3 2 )a^{k)a{k) + {a 2 - p 2 )b\k)b{k)) 
2n 

- (3 2 )(-N + N). (27) 



The energy momentum tensor 

Tyv = tpil^ip - T] ftv ^ij a d a ^ (28) 

is likewise only non-zero for /x, v = 0, 3. The various components are related: Too = T33 = 
— nTo3 and lead to the zero mode contribution to the energy per unit length of the string 
(respectively to the tension): 

f°° dk 

Ep = —Tp = / — k(aUk)a(k) + 6 t (fc)6(fc)), (29) 
Jo 27r 

The quantisation of the Majorana field goes analogous. We write 

/°° r\k 
— [a{nk)^ n e- ik ^ + a\nk)^-^) , (30) 

supplemented by the anti-commutation relation {a(k), a'(k')} = 2ir8(k—k'). The creation 
and annihilation operators for the particle and the antiparticle are now identical. Formally 
the Majorana zero mode can be obtained from the Dirac zero mode by setting (3± — > =F«± 
and b(k) — > a(k). 

The calculation of the conserved currents gives the following. As for any Majorana 
particle, the vector current vanishes (a Majorana mass term is not invariant under a vector 
U(l)). The axial current is also zero, by virtue of Eq. ()23|) : the amount of left-chiral and 
right-chiral fields exited for each zero mode is equal. We thus find 

Q = 0, Q 5 = (31) 

This can also be found from Eq. (|27jl. by substituting b(k) — ► a(k) and (3 — > a, as 
appropriate for a Majorana spinor. However, the energy-momentum tensor is non-zero 
for /i, v = 0, 3, as can be seen from Eqs. (|2*Hl [2HJ). This gives a non-zero contribution to 
the string energy per-unit-length and the string tension proportional to the number of 
Majorana particles trapped on the string: 

r au 

Ep = —Tp = / — k(a{k)^a{k)). (32) 
Jk>o 27r 

Before we go to a discussion of the effective 2 dimensional theory, a word on nota- 
tion/nomenclature. In the following we will refer to Q as particle number, the number 
of particles minus anti-particles (it is zero in the Majorana case). The term excitation 
number will be reserved for the total number of excitations, particles plus anti-particles. 
The excitation number does not correspond to a conserved current. It can be inferred 
from energy-momentum. 
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3.3 Dimensional reduction 



The zero-mode solutions effectively live in 2 dimensions. In this section we discuss how 
to dimensionally reduce from 4 to 2 dimensions. 
Define the matrices 

r 0± = ^(±7° + 7 3 ), r 1± = i(7 l ±fy 2 ). (33) 

They satisfy the algebra {P + ,r j ~} = 5^ and {P^r- 71 * 1 } = 0, and thus act as creation 
and annihilation operators. Define the ground state ( as the state which is annihilated by 
= for all i. We can then construct all states in the spinor representation by acting 
with creation operators on the vacuum. In 4D these states can be labelled by (s ,Si), 
with Si = ±1/2, corresponding to the state 

= /pO+-jso+l/2^pl+-jsi+l/2^ (34) 

The zero mode solutions , see Eq. (flHjl. in this language are 

< = a+(-l/2,-l/2)+/3+(-l/2,+l/2) 

tfj°_ = a_(+l/2,+l/2) + /3_(+l/2,-l/2). (35) 

Under dimensional reduction from 4 to 2 dimensions 

50(3, 1) -> 50(1, 1) x SO (2) ~ 50(1, 1) x U(l) (36) 

the spinor decomposes as 4 — ► 2 + 2. The 2's are spinor representation of 50(1, 1), which 
are charged under the action of the global U(l). We can construct the 2D spinor represen- 
tation analogous to the 4D construction. There is only one pair of creation/annihilation 
operators, given by 

f 0± = i(±7° + 7 1 ). (37) 

Here and in the following we denote all 2D quantities by a tilde. For the 2D gamma 
matrices we use 7 = o\ and 7 1 = ia-x- The 2D spinors are labelled by (sq). Explic- 
itly, (!) = (1 0) T which corresponds to a right moving fermion, and (— |) = (0 1) T 
corresponding to a left moving fermion. 

We dimensionally reduce (sq, s\) — > (sq) Si , so that right/left moving in 4D corresponds 
to right/left moving in 2D. The s\ labels the charge under the global U(l). 3 Applying 
3 This can formally be implemented by introducing a matrix U — (I2 o"i) T such that U'T Q+ U = T 0+ . 
Then tp = Uip is the same operation as (so, si) — > (so)- Helicity eigenstates, i.e., eigenstates of 7°7 3 7 5 i/' = 
±ijj, have si = ±1/2. 
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this to the zero mode wave function we get 

Each zero mode has 2 degrees of freedom (particle/antiparticle), and can be described by 
a single 2D complex Weyl fermion. 
The quantised Weyl spinor is 

^ (iP° n ~a(nk)e- tk ( z - nt) + ^J(nk)e^) (39) 

with ip® given in Eq. (|38|). We can now compute the 2D conserved charges. It is important 
to distinguish between 2D and 4D chirality. In D dimensions, chirality refers to left and 
right-handed particles, eigenstates of Pl,r = (1d±7d)/2. The zero mode solution is 
a superposition of both 4D chiralities. In 2D 7 s = diag(— 1,1). Hence, 2D chirality is 
equivalent to left and right moving particles, which is unrelated to the chirality concept 
in 4D. In 2D the vector and axial currents are related j = e^j". The two conserved 
charges are 

Q n=1 = N L -N L and Q 5 n=1 = -{N L - N L ) 
or Q n= _ l = N k -N k and Ql = _ 1 = N k -N R 

with L, R denoting 2D chiralities, i.e., left/right moving. In 4D left/right moving is related 
to fermion getting their mass via coupling to <fi or 0*. 2D particle number Q corresponds 
to 4D particle number Q. The 4D axial current corresponds to the 2D global U(l) current, 
under which states with different si are charged differently. Furthermore, the 2D energy 
and momentum are equivalent to the 4D expressions. 

The 2D Majorana fermion is obtained by setting a = (3 in Eq. ([38)1 . It has 1 real 
d.o.f., corresponding to a Majorana- Weyl fermion. The 2D charge conjugation matrix is 
C = I2. In addition, under charge conjugation also si «-> — s\. There is no global U(l) 
symmetry, no vector or axial charges but there is a non-zero energy and momentum. 

3.4 Bosonisation 

If is well know that a 2D fermionic theory with Dirac fermions is equivalent to a 2D bosonic 
theory. Here, we remind the reader how fermionic superconductivity can be discussed in 
terms of bosonic superconductivity ^T] . 
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(38) 



If there are both left and right moving fermionic zero modes the theory is anomaly 
free if the sum of charges of the left moving particles is equal to those of the right mov- 
ing particles. A theory with two zero modes of 4D Dirac particles moving in opposite 
directions is effectively equivalent to two 2D Weyl fermions with opposite chirality with 
the same charge. These combine to form a single 2D Dirac fermion. Such a theory can 
be bosonised. If the fermion is charged with charge e, it is equivalent to a massive scalar 
boson with mass e/y/n. 

If on the other hand there is only a left (or right) moving Dirac zero mode, which is 
necessarily neutral, the effective theory consists of a single 2D Weyl fermion which has 2 
real degrees of freedom. It is equivalent to a free massless chiral boson with an action 

C B = \d^dy. (41) 

subject to the constraint: 

(d t ± d z )<j> = 0. (42) 

The bosonic action in Eq. ()41)1 is invariant under <fi — > <fi + n, with n integer (the 
integer comes from the integer valued fermion current). Hence, is an angular variable. 
A topological invariant can be constructed for a string loop Q = J dl(d z <j)) [TT]. The 
topological invariant Q corresponds to particle number Q = N — N of the fermionic 
system. 

A Majorana zero mode has only one real d.o.f and can be described by a 2D Majorana- 
Weyl fermion. Bosonisation is not possible, as a Major ana- Weyl spinor has central charge 
c = 1/2, whereas a chiral boson has central charge c = 1 — the degrees of freedom do not 
match. Vortons studies are commonly done in terms of the bosonic description and it is 
not clear that this is adequate for a Majorana. 

3.5 Discussion 

The difference between the neutral Dirac current and the Majorana current is the follow- 
ing. 

For the Dirac case particle number is non zero and flowing along the string. It is a 
conserved quantum number. Current conservation is guaranteed by an (accidental) U(l) 
quantum number (such as lepton number). Two leptons cannot pair annihilate, and one 
lepton cannot decay into a state without lepton number. A left (right) moving zero mode 
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can be described by a left (right) moving Weyl fermion. The action can be bosonised, 
and is equivalent to a chiral boson. One can define a winding number which is conserved. 
This topological invariant can be defined for the boson field which corresponds to particle 
number in the fermionic description. 

If there is only one fermionic zero mode the string is chiral. This can only happen 
for neutral fermions, otherwise there would be an anomaly. Although the current might 
be persistent, it does not correspond to superconductivity, in the sense that there is no 
Meisner effect. Moreover, in the fermionic language, there are no cooper pairs. Cooper 
pairs are bound states of fermions with opposite spin and momenta, which therefore can 
only form if there are both left and right moving currents [TH] . 

For the Majorana case, there is no chiral or vector charge. Thus no-topological invari- 
ant can be defined. This is related to the fact that a 2D theory with an odd number of 
Majorana- Weyl fermions cannot be bosonised. Particle number is zero but in principle 
there can be a non-vanishing number of particles: there is then an energy-momentum 
flow. Since the Majorana does not carry a conserved quantum number, the number of 
Majorana excitations is not a conserved quantity in general. However, massless excita- 
tions cannot decay, neither pair annihilate (as they do not scatter), and thus the current 
is stable for kinematical reasons in the massless limit. 

Note that the masslessness of the zero mode is only a good approximation for an 
infinitely long straight string. String curvature and string oscillations will induce an 
effective mass. This will be discussed in more detail in the next section in the context 
of vorton formation. The other issues we want to address is whether a non negligible 
amount of chiral fermions can actually be trapped as transverse zero modes on the string, 
and contribute to the string energy-tension, and what is the actual difference between a 
flow of Majorana- Weyl fermions and a current of neutral Weyl fermions. 

4 Vortons 

In this section we discuss fermionic superconductivity in cosmic strings. The first ingre- 
dient we need is a fermionic zero mode. In the previous section, we have discussed the 
zero mode in detail for both Dirac and Majorana particles. Once we have a model which 
admits zero mode solutions, we need to know whether these zero modes can be exited or 
can be trapped so as to generate a current. Then finally, we must check for the classical 
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and quantum stabilities of these currents. 

In the case of charged particles, a current can be generated as the string moves through 
a magnetic field [H]. This is not possible for neutral particles. The only way that a large 
loop current can be generated is if there is an initial tiny current. As the loop collapses, 
the line integral of the current remains constant so that the current grows as the inverse 
of the loop length. In the case of charged particles, there are both left and right moving 
modes in order to cancel anomaly; charged fermions can thus scatter and this process 
strongly limits the current [13 j • In the case of a neutral particle, there is only left or right 
moving modes, and the scattering process is absent |13j . Therefore large currents may 
only occur if there is an asymmetry between left and right moving modes. This chiral 
asymmetry can only happen for neutral fermions which can either be Dirac or Majorana. 

Neutral currents cannot induce any (large scale) electric or magnetic field. Their 
cosmological interest lies solely in the conjecture that fermionic currents can stabilise 
string loops against gravitational collapse. Such stable, superconducting string loops 
are called vortons ^2]- Vortons contribute to the dark matter in the universe. The 
requirement that they do not overdose the universe can put strong constraints on the 
underlying particle physics [T§] . 

Vortons can be understood semi-classically [T2] HI] . Consider a fermionic current 
flowing along a string loop. Assume that the loop curvature does not significantly alter 
the zero mode spectrum obtained in the limit of zero curvature. To simplify all formulas, 
we further assume that for the Dirac current iV ^> A. 4 The zero mode energy levels are 
e« = n/R, with R the radius of the string loop. Each level can be occupied by at most 
one fermion and one anti-fermion, in accordance with Pauli exclusion principle. In the 
ground state, the highest occupied level has Fermi energy e = N/R. The total fermionic 
energy is obtained by summing over all occupied levels. The total energy of string is the 
sum of energy in bosonic and fermionic excitations E = Eb + Ep: 

N 2 

E = 2nfxR+— (43) 

where we have set N(N + 1) — > iV 2 , valid in the limit N ^> 1. Further \i ~ 0q is the 

string energy per unit length. As the loops shrinks under gravitational action the bosonic 

4 For N ~ N we find that the Fermi level is e = max[N,N]/R and fermionic energy Ej — (N 2 + 
N 2 )/(2R), results which only differ by order unity from those for N ^> N. In the literature vortons are 
often discussed in bosonised language, in which the only quantum number is particle number Q = N — N, 
which can be substantially different from N + N. 
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energy decreases, but the fermionic energy increases due to an increasing spacing between 
energy levels. The vorton radius, the radius for which the energy is minimised, is 

* = vfc (44) 

The above semi-classical description is valid for a loop radius larger than the string width 
and than the Compton wavelength of the zero mode. In particular, the last requirement, 
R Q > m^ 1 , gives 

N > A -1 , (45) 

with A the Yukawa coupling, see Eqs. (fHfl I2U|) . We will discuss current formation and 
stability for Dirac and Majorana currents, which may result in vorton formation, in the 
rest of this section. 



4.1 Effective zero mode mass 

To understand vorton formation, to understand current formation and stability, it is 
important to realise that the strings in the early universe are not infinitely long straight 
strings. Rather, the strings vibrate and curve, there are cusps and intersections. The 
zero-mode, which is a strictly massless excitation (in the sense that it travels at the speed 
of light) only for an infinitely long straight string, obtains an effective mass under such 
condition. Process which are absent for massless zero modes, decay is an obvious example, 
become possible. It goes without saying therefore, that an effective mass for the zero 
mode can alter the formation and stability of the fermionic currents considerably. Before 
considering these processes in more detail, we will first discuss under which conditions 
the zero mode gets effectively massive. 

The zero- mode lives in a potential well of height E = m^. At the bottom of the well 
the excitation is massless. What keeps the zero mode in a circular orbit for a string loop 
is that as the mode strays off the centre of the vortex (bottom of the well) it becomes 
massive. There is a restoring force (the centripetal force) —dV/dr ~ —5m/5r to keep it 
in orbit. The effective mass of the zero mode excitation can be estimated as follows [2U] . 
The acceleration of a particle going at the speed of light around a circle with radius R 
is \k\ = k/R. Setting this equal to the centripetal force, and reminding that Sr ~ m^ 1 
with rriij} the vacuum fermion mass (since the zero mode wave function falls off as e~ m ^ r 
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at large r) we get 

k k / a \ 

m ~ ~ . (46) 

The higher the momentum of the mode the larger the effective mass. 

Another way to see that string curvature induces an effective mass (alters the disper- 
sion relation to E ^ \k\) is to consider the solutions to the fermionic equations of motion. 
Consider a string loop with radius R lying in the xy-plane, with its centre at the origin. 
Define also spherical coordinates (r 1 , 9') which parameterise the plane transverse to the 
string, with the origin at the string core. This plane is characterised by r' x 9' || k, and 
A±_ = 0. The zero mode solution is schematically of the form 

$ = e*V st+t *)il) (T>,&). (47) 

Note that both k and r' x 9' are functions of time, since as the zero mode travels along 
the string loop k/\k\ changes. As a result the dt operator in the equation of motion not 
only brings down a factor E, but also terms proportional to k/k. This gives a dispersion 
relation of the form 

E ~ k + \ (48) 

K 

Expanding E = \/k 2 + m 2 « k + m 2 /k valid for m <C k, and setting k ~ mm^ equal to 
the centripetal force, this gives the same mass correction as the classical argument leading 
to Eq. flUl). 

Another string configuration which leads to an effective mass for the zero modes is a 
vibrating string. Consider a string extended along the ^-direction, which moves with some 
velocity in the x-direction. The dispersion relation for the zero mode is E 2 = k^. + k 2 z . By 
performing a Lorentz transformation we can go to a frame in which the string is at rest, 
and E = k z . However, for a vibrating string we cannot go to a frame in which the string is 
at rest at all times, as k ^ 0. Just as for the case of a string loop, the dispersion relation 
picks up a term k/k, see Eq. (|48|). For string loops the natural vibration frequency is 
uj ~ 1/R, whereas for string vibrations the natural scale is to ~ yjji ~ <f)Q. The former 
case gives an effective mass as in Eq. (|46|) . whereas the latter gives a mass which is larger 
by a factor 4>qR. 

Finally, effective mass terms are generated if locally components of the gauge field 
other than Aq are excited, or if the Higgs field profile deviates from its cylindric symmetric 
form. The string profile can be disturbed in string intersections or in cusps. A non-zero 
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A z or A component leads to a dispersion relation of the form 

E = k + qA, (49) 

and the effective mass is m 2 = qAk. A non-zero A r (or non-cylindrical component of 
<p) will appear as an effective mass term in the r, 9 dependent part of the equation of 
motion. It can be removed from the equation of motion by multiplying the zero mode 
wave function by a factor e iqAz . It thus likewise gives a dispersion relation of the form 
Eq. (|49|). In cusps and intersections A~l/r~mjin Eq. ([49)1 . and the effective mass is 
large m ~ k. 

The fermionic currents back react on the gauge field. This creates an effective mass 
for non-chiral strings but not for chiral strings. 

4.2 Current formation 

We consider chiral strings only (as appropriate for Majorana currents) and thus chargeless 
current carriers only. There is no possibility of current creation through spectral flow, as 
in the original Witten model, as the fermionic zero modes are neutral. Hence, any net 
current must be the result of fermions getting trapped/condensing on the string. Currents 
may be created in thermal or non-thermal settings. 

Thermal trapping. Let's first consider the case in which the string is emerged in a 
thermal bath with a reheat temperature larger than the vacuum fermion mass Tr > m^. 
Whatever the specifics of the trapping mechanism, it is expected to be inefficient at high 
temperature T 3> m^. The reason is that the typical energies involved in the trapping 
reactions are much higher than the binding energy of the zero mode Ef, = —m^, and 
trapping is suppressed. Fermions that do get captured by the string have typical momenta 
k ^> m^p and can escape easily at cusps and intersections. Furthermore, any zero mode 
present with k < can be easily scattered off the string. Hence, at high temperature 
we expect the number of zero modes to be negligible small. Current condensation has to 
take place at temperatures T < m^. 

On the other hand, at temperatures T <C capturing processes are also suppressed. 
For processes with one or more free fermions ip m t ne in-state this is obvious, as their 
number density is Boltzmann suppressed. But also if the initial state contains only light 
particles there is a suppressions factor. The reason is that the typical wavelength probed 
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in the reaction is much larger than the width of the zero mode, and the overlap between the 
zero mode wave function and those of the initial particles is small. As a result reaction 
rates are smaller by a factor (E/m^) 2 compared to the rate expected on dimensional 
grounds, with E ~ T the typical interaction energy (see appendix [2TJ 122 ) ■ 

We thus conclude that current formation, interactions in which zero modes are cap- 
tured/produced, has to take place at ~ T. 

Let us now consider the possible processes in which zero modes can be produced. The 
first interaction that comes to mind is a fermion scattering off the string to produce a 
zero mode final state, that is, the reaction ip — ► ip zm in the string background. The string 
is just a static background, it cannot exchange longitudinal momentum or spin with the 
fermionic states. Therefore energy and longitudinal momentum cannot be conserved in 
the reaction if the mass of the incoming fermion is larger than the effective mass m of the 
zero mode, see Eq. (jl(J8|) . And thus zero mode production in fermion-string scattering 
does not occur. 

But there are other reactions possible through which zero modes can be produced. If 
the only interactions present in the theory are those given in Eqs. (fEfl I20J1 . these are of 
the form ip + ip — > tp zm + if). 5 

In the limit T ~ the capturing and scattering rate are of the same order of 
magnitude, and if the reaction is sufficiently fast, a thermal equilibrium distribution of 
zero modes will be established. We thus define the number of zero modes on a loop of 
length L by 

N~N + N = P j. (50) 

with £ ~ T _1 the thermal correlation length. The efficiency factor p is of order unity if 
thermal equilibrium is reached, and p ~ if trapping reactions are frozen out at T ~ m/. 
The average particle number over many loops is zero (N — N) = 0. However, the variance, 
which gives the order of magnitude of particle number on the individual loop is non-zero: 

a/ ((N — N) 2 ) ~ y/N. 

Thermal equilibrium can be established if the reaction rate is sufficiently fast compared 
to the Hubble rate T ~ n^a > H . The thermal fermion density is ~ T 3 ~ m% at 
the temperature of interest. The cross section for the four-fermion interaction at energies 

5 An interaction such as tp + <fi — > ip zm , possible at T ~ if < m~0 has infinitely small phase space 
as the spin of the incoming fermion has to match that of the zero mode. 
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T ~ m^ is derived in appendix ID. 21 It is 

a ~ g * x A 6 ^ -. (51) 

max[mx, Wty] 

As discussed in the appendix we expect the amplification factor A to be of order unity. 
If the reaction is mediated by the Higgs field m x = m^ and gx = A is the Yukawa, 
whereas if the reaction is gauge mediated we have mx = mA and g x = g the fermion 
gauge coupling. The reaction rate exceeds the Hubble rate for 

3 x 10 16 GeV for m^ ~ 10 16 GeV 
"' V < {mi Mpi) 1 / 4 = \ 1 x 10 15 GeV for m^ ~ 10 14 GeV (52) 



g x AW 



3 x 10 13 GeV for m^ ~ 10 12 GeV 



If mx < m^, one should replace m x — > m^ in the above equation. If indeed A ~ 1 and 
gx ^$ 1 this becomes impossible as m^ — > 10 16 GeV. 

It is possible that the zero mode fermion ip couples to lighter fermions \i f° r example, 
through an interaction Cj = hHip^XL + h.c. This is of the form of the coupling of 
the right-handed neutrino to the left-handed one. Such a coupling allows for trapping 
reactions \ + X ~ > + 4> mediated by a light Higgs field. 6 At energies T ~ m^ 
(threshold energy) the cross section is a ~ h 4 mj 2 . This is of course of the same form as 
the cross section Eq. (joTj) in the limit mx < mj, and thus the bound on h derived below 
applies also to gx in this limit. The trapping rate is T ~ h^m^, which exceeds the Hubble 
rate at temperatures T ~ m^ for 



h > [jf-] = < 0.1 for m^ ~ 10 14 GeV (53) 
' 10 12 GeV 

Large couplings are needed for effective current formation. 

We conclude that thermal current condensation can only takes place at T ~ m^, and 

only if the trapping reactions are fast enough. This last requirement is hard to satisfy in 

the limit m^ < mA, m^ and in the absence of large Yukawa couplings to light fermions. If 

large Yukawa are present in the theory, current formation is expected to happen for m^ 

sufficiently small, with the number of zero modes given by Eq. (|5()jl . 

6 The reaction /' + h —> f zm is phase space suppressed, as the incoming light fermion needs to have 
the same spin as the zero mode. 
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1/4 




10 16 GeV 



For Dirac spinors lepton number is substantially less than the the total number of 
zero modes: Q = N — N ~ y/N + N. Since massless particles do not annihilate, the 
current description should be in terms of excitation rather than particle number. In the 
literature, vorton formation is often discussed in bosonised language in which the only 
quantum number is lepton number. The bosonic language strictly only applies if all 
particle anti-particle pairs have already annihilated, which would mean zero current for 
the Majorana case. 

Non-thermal capturing Fermionic current condensation can also happen non-thermally. 
If there is a large amount of energy released inside the string, this will create bosonic exci- 
tations (A and 0), which decay into fermions. If decay is immediate, T^c > H evaluated at 
the time of the energy release, the bosons decay inside the string, and the decay fermions 
may be captured as zero-modes. This was first discussed in [21], where they assume the 
energy-release is due to an internal (only inside the string) phase transition. This needs 
a complicated Higgs sector, and we will not discuss it further. 

One could also consider the case T R < m^. Before decay the Higgs field oscillates in 
its potential. If w < w 2 with w = k 2 + m 2 s ((f)) for all particles coupling to the Higgs field 
(0, A, x/j), then the oscillation is adiabatic. The string profile function and the fermion 
zero modes change adiabatically to adapt to the changing background fields. However, 
in the opposite limit, the evolution is non-adiabatically. Parametric resonance leads to 
explosive particle production. The fermions thus produced inside the string (from direct 
amplification, or from decay of the bosonic fields) can be trapped. It is not clear whether 
this mechanism can be efficient. One suppression factor is that only soft quanta are ex- 
ited during preheating with typical wave numbers much larger than the string width. 
Trapping is then suppressed by factors (k^/m^). 

To summarise, chiral currents which are necessarily neutral cannot be created by 
spectral flow. Current formation can still take place, either thermally or non-thermally. 
However, both processes are not obvious, and require special circumstances. 

4.3 Current stability 

The effect of an (momentum dependent) effective mass opens up the possibility for zero 
modes leaving the string, for pair annihilation, and for zero mode decay. We will discuss 
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these processes in turn. 

Classically, the zero mode cannot leave the string as long as m < m^. Conservation 
of angular momentum requires the massive mode outside to take away all the momentum 
of the zero mode. Consequently, there is an energy barrier if m < m^, but not in the 
opposite limit. As first noted in [20], the generation of an effective mass as in Eq. ()46j) 
puts a maximum on the current. For a vorton loop m ~ k/(Ram^) with Rq given in 
Eq. (pFfj) . Using that the zero mode momentum is typically of the order of the Fermi 
momentum kf ~ N/Rq, the typical mass is 

m ~**&~ 10-GeV (^) M^Y f^S^L) . (54) 
Requiring m < so that the current is classically stable gives 

\ 2 N > 4tt, (55) 

where we have used that = A0o- For perturbative couplings this gives a stronger 
constraint than Eq. (|45j). Although classically the zero mode cannot leave the string 
if m < rrty, there is a non-zero tunnelling probability for this process to happen [T4*] . 
The extra energy needed to overcome the energy barrier can be provided by the string: 
enlarging the loop radius will increase the bosonic energy, see Eq. (|4Hj) . The current is 
stable against tunnelling processes (lifetime longer than age of the universe) for 



with a = g 2 /(47r) the fine-structure constant of the broken U(l)x- 

Another consequence of a non-zero mass is the possibility of zero mode decay. As 
we have seen, thermal current formation is most effective if the zero mode has a large, 
order unity, coupling to light particles. Decay is kinematically allowed if the zero mode 
is heavier than the decay products. The decay rate for a coupling h, see Eq 1)990. then is 




r ~ -^2- I 1 + °\ — ) I 

For m the Compton wavelength m _1 is much smaller than the width of the zero 

mode, and the wave function overlap of the zero mode and the decay products is small. 
This gives a suppression factor (m/m^) 2 . For order one couplings and masses m > 10 GeV, 
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decay occurs before big bang nucleosynthesis (BBN). It follows that current stability can 
only be assured if the decay is kinematically forbidden and m < ma ec . 
A lower bound on the decay width comes from gravitational decay 



where we have again included a suppression factor (m/m^) 2 . Decay happens before BBN, 
and thereby evading all bounds from vorton overdensities, if 



This is generically satisfied for high scale phase transitions unless N is extremely large. 
Decay of Majorana particles cannot be suppressed by invoking symmetries; decay of Dirac 
particles can be suppressed by lepton number conservation, although it is believed that 
gravity violates all global symmetries. 

A third effect of a non-zero mass is that different modes travel at different velocities and 
can scatter off each other. Particles and anti-particle states (Dirac case) or two particles 
can pair annihilate (Majorana case) into decay products whose masses are lighter than 
the Fermi energy. The typical centre of mass energy is of the order of the zero mode mass. 
Thus annihilation is only possible if the mass of one of the zero modes is larger than the 
masses of the decay products. If Eq. (}55|) is violated, not only can zero modes leave the 
string, also the reaction ip zm + ip zm — > ip + ip becomes possible. Scattering into lighter 
particles, ip zm + ip zm — > X is kinematically accessible for nix < m^. The conditions under 
which scattering is effective are similar to those under which zero mode decay occurs. 
Note that pair annihilation into light particles cannot be forbidden by lepton number. If 
it does happen, the Majorana current is destroyed, whereas a Dirac current can subsist if 
there is a net lepton number. 

To summarise, the zero-modes obtain an effective mass from string loop curvature, 
string vibrations, and string perturbations making escape, pair annihilation and decay of 
the zero modes possible. Escape can be prohibited for large enough N. Annihilation and 
decay processes might be suppressed or forbidden for Dirac fermions by conservation of 
lepton number (though gravity is believed to violate global quantum numbers), but not 
for Majorana fermions. Decay and annihilation can occur at relatively low Fermi levels 
if the zero modes couple to light fermions — which is precisely the favoured condition 




(58) 




(59) 
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for thermal current formation. But even in the absence of light fermions, gravitational 
decay happens before nucleosynthesis unless N is very large or the string forming phase 
transition happens at low scale. Vorton formation does not look probable, especially for 
Majorana fermions. 

5 Conclusions 

In this paper we have derived the zero mode solution for a Majorana neutrino and com- 
pared it with its Dirac counterpart. The Majorana zero mode has one degree of freedom, 
one less than the Dirac zero mode. Both the vector and axial currents are zero. It is 
therefore impossible to bosonise the 2D effective action for a Majorana zero mode. 

No conserved quantum number can be defined to assure the stability of a Majorana zero 
mode. This is in contrast with Dirac fermions, for which always the analogue of lepton 
number can be introduced. However, massless zero modes are stable for kinematical 
reasons. And in the massless limits there is thus no difference between the stability of 
Majorana and Dirac zero modes. 

The zero mode is only strictly massless for infinitely long straight strings. The cosmic 
strings in our universe are expected to be curved, to vibrate, to fold and to intersect. 
Under such conditions an effective mass for the zero mode is generated. Needless to say 
that an effective mass alters the stability properties of a zero mode. We have studied 
these issues in the context of vortons, string loops stabilised by a fermionic current. 

Neutral fermionic currents can only be formed if zero modes are somehow trapped on 
the string. A massive fermion outside cannot simply be trapped through scattering with 
the string, as this would violate energy-momentum conservation. Interactions in which 
more than two initial/final state particles participate are needed. Such interactions can 
only be efficient if the zero mode couples to light particles with large, order one, couplings. 

The fermionic current on a string loop increases as the loop contracts. The effective 
zero mode mass, which is proportional to the loop curvature, increases in this process, 
opening up the possibility of zero mode decay. The decay width is large if the zero 
mode has a large coupling to light particles, precisely the conditions under which current 
formation is effective. We expect that vortons will not form unless the string scale is low. 
This is especially true for Majorana neutrinos which can pair annihilate, and whose decay 
cannot be protected by a conserved quantum number. 
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A Notation 



We use the Minkowski metric rj^ = diag(l, — 1, —1, —1). For the gamma matrices we take 
chiral basis: 

o)' " = {o J (60) 

with <7 M = (1, cr l ) and <x M = (1, — a 1 ), and a 1 the Pauli matrices: 

•*-(:;)■ -(::•) <«> 

In this basis, the projection operators Pl,r are defined by 

^ = ^L+^fl = ^ + ^*= ^T~^ + ^T~^ ( 62 ) 

The charge conjugate of a spinor is defined as 

^ = cy,t = and ^ = ^ 

with 

C = n\ 

C = C( 7 °) r = *7 2 7°- (64) 

The latter matrix has the properties C = —C^ 1 = —C^ = —C T . The charge conjugate 
field ifj c acts similarly as the hermitian conjugate field = ip^'j : it creates particles and 
annihilates anti-particles. The chiral projections of ip c are given by 

r L ,R = C$l L and 4>l R = V R , L C (65) 

Note that ip c L = (ip c ) L = P L *p c = P L C$ T = CP L $ T = C(ipP L ) T = C$1 = (ip R ) c is a 
left-handed spinor, not to be confused with {jPl) c = {Pl$)° = CVl = (i/> c )r = V'a "which 
is a right-handed spinor. 



B Spinor conventions 



In this appendix we list our spinor conventions. The normalization is slightly different 
from that used in the paper. 
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B.l Free Dirac spinor 

We first give the results for a free massive Dirac spinor. The positive and negative 
frequency solutions of the Dirac equation are 

with r = ± corresponding the two independent helicity/spin states. In our conventions 
£ + = ( Q ) and = (°) . The spinors are normalized to 

u\p,s)u(p,r) = v^p, s)v(p,r) = 2E p S rs , (67) 

u(p, s)u(p,r) = —v(p,s)v(p,r) = 2mS rs , (68) 

u\p, s)v (— p, r) = v\—p ) s)u{p,r) =0, (69) 

u(p, s)v(p, r) = v(p,r)u(p,r) =0. (70) 



The quantized field is 



I Wf^kw ^ ( as{p)u{p > s)e ~ w ' x + b ^ p)v ^ s)eipx ^ ■ (71) 

Imposing an equal-time anti-commutation relation for tp gives the anti-commutation re- 
lations for the creation and annihilation operators: 

{^ f } = 5 3 (x-y)U {a s (p),at(k)} = {b s (p),bt(k)} = (2n) 3 5 3 (p- k)5 rs (72) 

with 1 4 the unity matrix. We define the one-particle state as 

H = y/2E p a\{p)\Q). (73) 

so that 

(ps\kr) = 2E p (2ir) 3 5 3 (p- k)5 rs (74) 

where we have used that the vacuum is normalized to unity (0 1 0) = 1. Our renormalization 
choice Eq. (JZHJ) corresponds to a wave function normalization of 2E P particles per unit 
volume. The wave function is 

<0|V|ps> = u{p,s)e- ipx (75) 
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B.2 Zero mode spinors 

The classical zero mode solution to the wave equation is (we consider only n = 1, and 
coupling to 0*) 

*>(!>)= frS) (76) 



This is the solution for both positive and negative frequencies. Further we require 



J dxdy (a 2 + (3 2 ) = 1 (77) 



so that the spinors are normalized 



dxdy V\pW\p) = 2E P 
tp°(p)ip°(p) = ip°\p)ip°(-p) = (78) 



We approximate the profile functions by 

M 



a = -(3 = ^e~ Mr , (79) 



7T 



satisfying the normalization condition Eq. (fTTj) . The zero modes for a vortex/anti- vortex 
with \n\ = 1 are independent of the angular variable 9. The width of the wave function is 
set by the fermion mass outside the string M = m^ 7 , as the solution to the Dirac equation 
at large r is oc e~ m ^ r . 

The quantized zero mode field is 

^ = j ^-J= M P V^-^ + «e'^) (80) 

with p = \pl\ = E p with pi the momentum longitudinal to the string. The field anti- 
commutator reads 

J dxdy =^ Hp),J(k)} = {6(p),6t(fc)} = (2n)5 1 L (p- k) (81) 

with 1 4 = diag(0, 1, 1,0). The delta function involves only the momentum aligned with 
the string. The matrix I4 instead of the unity matrix appears because the zero mode has 
only one spin state, there is no summation over two independent spin states as for the 



7 We introduce M so that in the interaction rate the origin of the different factors can be traced easily. 
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free Dirac field. The one-particle state following from it can be defined analogous to the 
free Dirac field Eqs. (J?3*| The one-particle state is 



\p,s) = ^2E p al(p)\0). 



(82) 



so that 



(ps\qr) = 2E p {2tt)5\p - k)5 r 



(83) 



The renormalization choice Eq. (|82J) corresponds to a wave function normalization of 2E P 
particles per unit length; it reads 



(0|VM = V{.e-** = u ^=e- Mr ^- x , 

\/vr 



(84) 



where we have defined u° through 



M 



-Mr 



71 



(85) 



The product u°u°, which is needed in the calculation of the squared amplitude \A\ 2 , in the 
frame in which the string is aligned with the z-axis and the zero mode has 4-momentum 
= (p,0,0,p), is 

/ \ 
-10 1 
10 0-1 
\ / 

Similar expressions exist for the anti-particle zero mode. 



u°u° = 2p 



(86) 



B.3 Interactions in the string background 

The interaction rate for some processes involving a zero mode will be calculated explicitly 
in the next two appendices. Here we will give some general remarks. 

Following fH\ E21 12H| interactions in the string background are calculated using the 
approximation that the incoming and outgoing non-bound states are asymptotically free. 
That is, we use a free plane wave expansion for these states. The plane wave expansion 
can then be matched to an angular mode expansion, which is more natural in the string 
background (as is done in Eq. (|9*K1 below) . As the zero mode is bound to the string, it is 
a good approximation to keep only the the lowest angular modes. Possible amplification 
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factors, due to an amplification of the fermion wave functions of the non-bound states 
near the string core, will be added by hand, see Eq. (jll4j) and further. 

The main difference between interaction rates involving only asymptotically free states, 
and interaction rates involving zero modes are the following. 

1. The zero mode has only one spin state. This gives factors of order one difference, 
but does not lead to large quantitative or qualitative differences. 

2. Conservation of transverse momentum is violated in the string background. For 
typical interaction energies E < this has no big consequences. 

3. The interaction rate is suppressed by a factor (E/M) 2 if the typical interaction 
energy E < M with M ~ the width of the zero mode. The reason is that the 
overlap between the zero mode wave function, and the wave functions of the other 
states is small [22]- 

4. Non bound states can have wave functions which are enhanced near the string core 
compared to the plane wave. This factor is model dependent as it depends on the 
core structure, the fractional flux and the fermion charges. This can lead to large 
amplifications of the reaction rate fU\ I23j. 

C Zero mode decay 

To set the ground and to introduce the notation we start with a discussion of the decay 
rate of a massless, asymptotically free fermion. 

C.l Decay of massless particle 

Consider the decay i/jn(k,r) — * H(p')xL,(k' ,r'), mediated by a term in the interaction 
Lagrangian 

Cj = hH^P LX + h.c. = hH$ RX L + h-c. (87) 

Decay is kinematically allowed only if the decay products are massless, and if their mo- 
menta are collinear (and aligned with the momentum of the decaying particle ipa). The 
phase space is infinitesimal small. Since all products p ■ k vanish, with p, k the 4-momenta 
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of the particles involved in the decay process, all Mandelstam variables are zero. The 
amplitude is 

A = J d 4 x (k'r';p'\£i\kr) 

= (-ih) J & A xu{k\r')P L u{k,r)e- i{k ~ k '- pl) 

= {2Txf5 i {k-k' -p')(-ih)u(k',r')P L u(k,r). (88) 

Here k is the momentum of the initial state fermion, k' that of the final state fermion, 
and p' that of the final state boson. We use box normalization. The interaction rate per 
unit volume can than be defined through 

p = \A\ 2 /VT = (2vr) 4 5 4 (A; — k' — p') \M\ 2 (89) 

with A = (27r) 4 <5 4 (fc — k! — p')iA4. The decay rate then is 

r = -j- / cWcLy(27r) 4 5 4 (fc - k' - p')\M\ 2 (90) 
2t, k J 

with dojp the invariant 3-volume 

dto p = - d " P (91) 

Averaging over initial spin states and summing over final states, the matrix element 
becomes 

h 2 

\M\ 2 = y y^Tr [Uk'r'PLUkrUkrPRUk'r'] 
r,r' 

= h 2 (k-k') 

= (92) 

where the last equality applies for massless collinear particles. Hence, the decay rate for 
a massless fermion is zero. 

C.2 Decay of the zero mode 

Consider now decay of the massless zero mode through the same channel Eq. (|87j). We 
have seen that the decay width for a massless free particle is proportional to k • k', which 
vanishes for massless collinear particles. Processes involving the zero mode do not conserve 
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momentum transverse to the string. This is a consequence of the fact that Lorentz 
symmetry (translations in the plane transverse to the string) is broken by the presence 
of the string. Therefore, the decay products of the zero modes can acquire transverse 
momentum and do not need to be collinear. Can this result in a non-zero decay width? 
Let's calculate. The matrix element for ip zm (k,r) — > tp(k',r') is 

{k'r'^PL^kr) = u k , r ,P L ^ e- i{k - k ' )LXL e- ik 'T XT 

= u k , s ,P L u° k e- i{k - k ' )LXL ^=e lk T XT - Mr (93) 



"K 

where in the last line we have used Eq. ()85j) . The subscript L and T denote space-time 
longitudinal (including the time-coordinate) and transverse to the string respectively. For 
the zero mode k T = 0. The decay amplitude for the process ip zm (k) — > H(p')ip(k', s') then 
is 



A = -ih J d 4 x(k'r'\i;P L i; zm \kr)(p'\H\0) 



= (2n) 2 5 2 (k - k' - p')(-ih)u k , r ,P L ul 9 -^ (94) 

with qjn = (k' + p')r- In the last line we have introduced the dimensionless function g 
through 

9(Qt) f j 2 „ M Ak> x T -Mr 



M J 

« ^2tt / drrJ (q T r)e- Mr = ^(1 + q %jM^ (95) 
V 71 " J M 

where in the first step we have only kept the lowest angular mode. There is no delta 
function imposing conservation of transverse momentum. Rather the integration over 
transverse coordinates lead to a distribution which has maximum g ~ 1\pn for qx = 
and decreases (power law) for 3> M. There is a non-zero probability for the decay 
products to acquire transverse momentum. 

Using box normalization we can define the interaction rate per unit length as 

P = \A\ 2 /LT = (2ti) 2 5 2 L (k - k' - p')\M\ 2 . (96) 

The decay width then is 

T = & k j d ^' d ^'^ 26 l( k ~ k '- P') (^fr) 2 E \ u ^Pl4\ 2 (97) 
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The delta function in the decay width can only be satisfied for massless decay products, 
with their longitudinal momentum in the same direction as the zero mode momentum. 
In particular, it implies Ey = k! and qx = 0. Since the decay product are then collinear 
just as in the decay of a non-bound fermion, we also expect the zero mode decay width 
to vanish. Indeed, the spinor factor, evaluated in the frame where the string is aligned 
with the z-axis so that the zero mode momentum is W = (k, 0, 0, k), is 

\uk>r>P L u° kr \ 2 = S s /Tr [u f P L u°u°P R u f ] = 2k{E k , - k' z ) = 0, (98) 

r' 

where in the second step we have used Eq. (JEHJ), and the equivalent for the massless free 
spinor = $ . 

Energy-momentum conservation in the longitudinal plane assures that a massless zero 
mode does not decay. 



C.3 Decay of massive bound state 

If the zero mode gets an effective mass m no longer E k = k, and thus no longer necessarily 
E' k = k', and the decay width can be non-zero. We will assume that the zero mode spinor 
is unaffected by the effective mass, only the frequency part is altered. The cross section 
is most easily evaluated in the center of mass frame, where the zero mode is at rest. 
The spin sum give J2 S ' I^v-Pl^J 2 ~ m2 - Introducing an extra delta function the decay 
width Eq. (|97|) then becomes 

„ h 2 m 



M 2 
h 2 m 



J dc^dovdV (27r) 2 c^(fc + q T -k' - p')g{q T f 

dV giqr) 2 



M 2 

h 2 m 3 / „/m\ 



O T7 ■ 09) 



M 2 \ \MJ 

Here, as before, qr = (k' + p')r is the transverse momentum provided by the string. The 
integration over k! and p' can be done in the usual way, and impose energy-momentum 
conservation in the longitudinal plane. Energy conservation restricts qx < m. In the last 
step we have taken the limit m«M~ m^, with the fermion vacuum mass, to arrive 
at the final result. For m -C M there is a wave function suppression. The reason is that 
then the Compton wavelength m _1 is much smaller than the width of the zero mode. 
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D Zero mode trapping 



We are interested in processes with in the initial state free fermion(s), and in the final 
state a zero mode. If such reactions are fast compared to the Hubble rate at T ~ m^, 
zero modes are trapped, and chiral currents form. 



D.l Zero mode from string scattering 

We first consider the process 

ip(k,r) — > ip zm (k',r') (100) 

in the background of the string. The interaction can take place via either the string 
Higgs or gauge field. As we will see, the reaction rate is zero as a consequence of energy- 
momentum conservation. 

The amplitude for ip(k, r) — > ip zm (k', r') through a fermion-Higgs coupling of the form 



C = X^iptp. (101) 

A = J d A x{S'-k'r'\d\S-kr) = J d 4 x(k'r'\^\kr)(S'\^\S) (102) 

with S denoting the string state. If the back-reaction is small, if the string is not much 

altered by the scattering, then S' S. Here $ = — O , the field shifted by the VEV 
outside the string. We approximate 

(S'\*\S) = I 00 ' r<m V (103) 
I 0, r > m , . 



is 



The amplitude then becomes 



A ~ (27i) 2 5 2 L (k - k')(-i\)u'°u5 rr/ ^ (104) 

m± 



where we have used that 



/ d 2 x T <j> -=e- ihl " Xt - Mr ~ M0 O / rdrJ (k T r)e 
Jo V n Jo 



Mr 



M*or i + (h.,*L }] (105) 



ml V V m </> m <t> 
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The integration over longitudinal momentum gives a delta function enforcing energy- 
momentum conservation in the longitudinal (t, z) plane. The string is just a stationary 
background, with no longitudinal string momentum. Note however, that the string can 
provide transverse momentum so that qx < E can be non-zero, with E the typical inter- 
action energy. 

Energy momentum conservation in the longitudinal plane, which reads explicitly 
Ek = E' k = k' and ks = k' 3 = k' , is only possible if the incoming fermion is also massless. 
The same conclusion follows from the Dirac delta in Eq. (jl04j) . The spin of the incoming 
fermion should be equal to the spin of the zero mode fermion. Only for massless incoming 
fermions is this possible. The massless spin states outside the string are chiral eigenstates 
and only in this case is it possible to write the zero mode wave function ^° as a superpo- 
sition of the fermion states outside the string. 8 We conclude that the process Eq. (jlOOj) 
to capture zero modes is impossible if the fermions are massive outside the string. 

We will show here that the cross section is also zero for massless incoming fermions. 
The cross section per unit length is 

da A 2 f dk' 1 - (MfoY nnR , 

di ~ E J (2^^ ~ k ¥-' |M M| [ijr) ■ (106) 

Evaluating the spin term in the frame where the zero mode momentum is k ,fl = (k', 0, 0, k') 
gives 

5 2 L {k - k')8 rr '\u°(k')u(k, s)\ 2 = 8\{k - k')2k'{E k - k 3 ) 

= (107) 

String scattering through the interaction term L = iqijj'y 6 Agip gives similar expression. 
It can easily be checked that also in this case the cross section is zero. 

Can string scattering lead to zero mode capture if the zero mode is effectively massive? 
Energy momentum conservation in the longitudinal plane gives k^+m^ = m 2 , with hp the 
transverse momentum of the incoming particle, the vacuum mass, and m the effective 
zero mode mass. It follows that zero mode capture through Eq. (jlOOj) is impossible for 

> m (108) 

which is the case of interest. 

8 The zero mode spinor is of the form u° = \/2k(0 1 1 0) T and the two massless free states are 
u + = V2k(0 1 0) T and w_ = V2k(0 1 0) T . It follows u° = u_ + u+. 
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D.2 Four fermion scattering 

Consider 

ip(k, r) + $(p, s) -> ^ zm {k') + tfj{p', s') (109) 

mediated by either a Higgs or gauge boson. Since the zero mode wave function For 
simplicity we consider only s-channel Higgs mediation explicitly; other diagrams and 
interference terms give similar contributions. The amplitude for this process is 

A = HA) 2 /'d 4 x 1 d 4 x 2 t;(p,g)M(fc,r) - 1 rf %k , )v(j/,3 , )e^ k ^ 9 > Xl ^ q -*- l/) - Xa - Mr 
J q 2 + mi 1 



(2n) 2 S 2 L (k + p-k>- p')^lv(p, s)u(k, r)u°(k')v(p', s' f {{q P ' )t) (110) 

ml M 



where we have assumed that the interaction takes place a distance r < mj 1 away from the 
string core. In the opposite limit r <C mj 1 the zero mode wave function approaches zero 
exponentially fast, and the cross section is zero. In the second line we have specialized to 
interaction energies q = k + p -C m^. The function g is defined in Eq. (195)1 . The cross 
section (for interactions taking place within a distance r < mj 1 from the string center, is 

& ~ ~~i f du p >du k ,d 2 q T Sl(q - k! - p')\T, rss/ \vuu°'v'\ 2 ^-^ (111) 
J 4 M z 

with day given in Eq. (|9T|) . and dujk' = dkj (47r/c). The spin sum is evaluated in the frame 
in which the zero-mode momentum is k ,fl = (k', 0, 0, k')\ 

^ r ^\vuu a 'v'\ 2 = Ak' ((k - k 3 )(m 2 -p-p') - CPq -Ps)(m 2 ~ k ■ p)) (112) 

which is non-zero for massive fermions (so that \k$\ ^ \k$\ and \p' Q \ ^ \p' 3 \)- For interaction 
energies q ~ the cross section then is 

4 f 2 
\4 m ^ / j2 2 M ~ m ^ \A m it> / 11Q \ 

a ~ X ^APj dqTg -> A ^ (113) 

For m^, -C M there is a (m^/M) 2 suppression as the Compton wavelength m^ 1 is much 
smaller than the width of the zero mode. In this calculation we have use that the free 
particles are of the plane wave form. However their wave function may be amplified on 
the scale of the typical interaction length m^ 1 . Defining 

\u(r = ml 1 )! 

^ = A ^ (H4) 

\u{r = m^)\ 
36 



with u the wave function in the string background and u, as before, the plane wave. The 
cross section corrected by these amplification factors becomes 

2 

a amp ~A 6 A 4 ^. (115) 

The amplification factor depends on the core profile and on the fractional flux. We expect 
the amplification factor to be 

A~(^Y<(^Ul (116) 



with the coefficient p < 1} 



kr I ~ V k 
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